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Abstract: We extend the range of possible applications of type analyses to decay 
chains with multiple invisible particles, as well as to asymmetric event topologies with 
different parent and/or different children particles. We advocate two possible approaches. 
In the first, we introduce suitably defined 3 + 1-dimensional analogues of the Mt2 variable, 
which take into account all relevant on-shell kinematic constraints in a given event topology. 
The second approach utilizes the conventional Mt2 variable, but its kinematic endpoint is 
suitably reinterpreted on a case by case basis, depending on the specific event topology at 
hand. We provide the general prescription for this reinterpretation, including the formulas 
relating the measured Mt2 endpoint (as a function of the test masses of all the invisible 
particles) to the underlying physical mass spectrum. We also provide analytical formulas 
for the shape of the differential distribution of the doubly projected Mt2j_ variable for 
the ten possible event topologies with one visible particle and up to two invisible particles 
per decay chain. We illustrate our results with the example of leptonic chargino decays 
X'^ — )• i^i^x^ in supersymmetry. 
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1 Introduction 

The Large Hadron Collider (LHC) at CERN is continuing its quest for new physics Beyond 
the Standard Model (BSM). Among the multitude of possible BSM scenarios, models 
with neutral and stable WIMPs (weakly interacting massive particles) are of particular 
interest. First, such models are greatly motivated by the dark matter problem, as WIMPs 
are suitable dark matter candidates. Second, the initial BSM searches at the LHC have 
already placed stringent limits on heavy resonances which decay visibly and can be fully 
reconstructed. In contrast, the limits on parent particles which decay semi-invisibly (to 
a collection of visible SM particles and one or more WIMPs) are much weaker. First, 
the background issue in these analyses is more complicated: since the parent cannot be 
fully reconstructed, the search is not a mere "bump hunt", where the background can be 
simply subtracted from the side-bands. Second, the symmetry which protects the lifetime 
of the WIMP dark matter candidate typically requires that the new particles are multiply 
produced, leading to lower production cross-sections (as opposed to single production). In 
the simplest and most popular models, the new symmetry is a Z2 parity, which implies 
that the new particles are pair produced, as illustrated in Fig. 1. 
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Figure 1. (a) The generic event topology for the pair production of two parent particles A and 
A' . Particle A {A') decays to a set of Nyis (iV^^) particles which are visible in the detector and 
a set of Ninv invisible particles. The invisible particles are not necessarily limited to SM 

neutrinos, and can have arbitrary masses, (b) The simplified event topology typically used in Mt2 
studies. Each V is either a single SM particle, or an effective composite visible particle which is 
constructed from the corresponding visible particles in figure (a). The decay chain for each parent 
A terminates in a single invisible particle C. 



Fig. 1(a) depicts the generic event topology in a typical SUSY-like missing energy event. 
Two parent particles A and A' are produced, and each one decays to a certain number of 
visible particles (denoted by solid lines), and a certain number of invisible particles (denoted 
by dashed lines). At this point, the diagram in Fig. 1(a) is meant to represent the most 
general case: for example, the number of visible particles {Nyis and -/V^g), as well as the 
number of invisible particles {Ninv and N-^^) in each parent decay chain is completely 
arbitrary. Furthermore, the invisible particles are not necessarily limited to SM neutrinos, 
but may include (several species of) WIMPs with different masses. Most importantly, 
no assumption has been made regarding the actual decay topology of the parents, and 
the yellow-shaded circles are simply placeholders for the actual 1 — )• Nms + Nnv and 
1 — )• N^ig + N^j^^ Feynman diagrams responsible for the parent decays in the figure. 

The kinematic analysis of events of this type is rather challenging. There is a lot 
of missing (or a priori unknown) information: the number of invisible particles Ninv and 
N-^^, the momenta (and masses) of the invisible particles, the total invariant mass and 
longitudinal momentum of the parent system {A, A'), and the event topologies hiding 
behind the yellow-shaded circles in Fig. 1(a). These problems have inspired a lot of previous 
work in the literature on new methods for measuring the masses of the parents and the 
invisible daughters under these circumstances (for a recent review, see [1]). Among the 
different options existing in the literature, the invariant mass variables appear to be both 
useful and theoretically motivated [2]. Being Lorentz-invariant, they are insensitive to the 
a priori unknown longitudinal boost of the {A, A') system. A special subset of invariant 
mass variables which also shares this property, is given by the transverse invariant mass 
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variables, among which the Oxbridge stransverse mass Mt2 [3, 4] is a well-known example 
directly applicable to Fig. 1(a). It has been featured in analyses performed by the Tevatron 
and LHC experimental collaborations CDF, DO, CMS and ATLAS [5-9]. 

In the original proposal [3] , the "Oxbridge" Mt2 variable was applied to the case of 
direct slepton production, corresponding to the simple event topology shown in Fig. 1(b), 
which is a special case of the more general Fig. 1(a). In Fig. 1(b), two identical parents 
(A) are produced, and each one subsequently decays via a two-body decay to a single 
visible SM particle V and a single invisible child particle C. In spite of its simplicity, 
the event topology of Fig. 1(b) covers a number of interesting physics cases, e.g. slepton 
production [3, 10], squark production [11, 12], chargino production [13, 14], etc. At the 
same time, the Mt2 concept is very powerful, and can be usefully applied to situations 
that are more general than the simple example in Fig. 1(b). 

So let us first review the different directions in which one could generalize Fig. 1(b). 
Some of the following options (a-b) have already been considered in the literature, and we 
only mention them here for completeness. The remaining possibilities (c-f), however, have 
attracted significantly less attention in the literature, and will be the main focus of this 
paper. 

(a) In Fig. 1(b) there is only one visible particle in each decay chain. In the language of 
Fig. 1(a), this implies the assumption Nyis = N^^g = 1. At the same type, a typical 
BSM model like supersymmetry exhibits much longer decay chains, with several 
visible particles on each side of the event. This case can be easily handled with the 
conventional Mt2 approach — one just needs to think of each ^ as a collection of 
visible particles with some net four- momentum P^, which is measured in the detector. 
Early work along these lines [15-17] led to the discovery of the Mt2 "kink" in the 
endpoint MJp^^ when considered as a function of the unknown test^ mass Mc for the 
child particle: 



The kink in eq. (1.1) is an interesting and unique property of the Mt2 variable, 
allowing to measure simultaneously (at least as a matter of principle) the true mass 
Mc of the invisible dark matter candidate C at the end of the decay chain and 
the true mass Ma of the parent particle initiating the decay chain. What makes 
the appearance of the kink possible is the fact that the invariant mass My of the 
visible collection of particles V is not constant, but varies from event to event: My = 
P^'P^ 7^ const. 

(b) Another possible generalization is to consider that the parent particles A are produced 
inclusively, either in association with jets from initial state radiation, or in the decays 
of other, even heavier, new particles. In either case, one again finds an Mt2 kink 
in eq. (1.1) at the true mass of the daughter particle C [17, 18]. The kink persists 

^Following the notation of [18], test input masses for invisible particles will be denoted by a tilde. 




(1.1) 
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even if the decay chain is extremely short {Nyis = N'^^^ = 1) and My = const. This 
is because the origin of the kink is now different — it is due to the net transverse 
momentum of the parent system {A, A). 

(c) Another hmiting assumption in Fig. 1(b) is that the two missing particles C at the 
end of each decay chain are the same (or at least have a common mass Mc)- This 
assumption, however, can be easily relaxed — one simply needs to allow for two 
independent mass inputs Mc^ and in the calculation of Mt2 [19]. The resulting 
"asymmetric" Mt2 variable inherits all the desired properties of the conventional 
Mt2- In particular, the "kink" in eq. (1.1) in the function MJp2^{Mc) is generalized 
to a kinky "crease" in the surface defined by the function MJp2^{Mci,Mc2) [19, 20]. 
While this all sounds very straightforward, one should keep in mind that the original 
public codes [21, 22] for calculating the Mt2 variable cannot be used in this case, 
since the assumption Mc^ = is already hardwired^. 

(d) The other assumption in Fig. 1(b) is that the two parent particles are the same (or at 
least have a common mass Ma)- In principle, this can also be handled rather easily. 
If the parents are different, but the children are the same, then one of the decay 
chains must have additional visible particles which are not present on the other side. 
Then, one possibility is to try to identify the extra particles and remove them from 
consideration, thus reducing the effective event topology back to Fig. 1(b) [24]. An 
alternative strategy is to consider the full event, but allow for different parent masses 
Mai and and then use M^^"" to construct the function Mc{Ma^,Ma2) [20]. 

(e) Another assumption of Fig. 1(b) is that there is a single invisible particle in each 
decay chain, i.e. that Nmv = ^Inv — 1- Again, there is no compelling reason for this 
restriction in light of generic BSM models. First, a Zi\i parity restricts the number 
of dark matter candidates in each decay chain only modulo N. Second, a Z3 parity 
can lead to Ninv = 2 as easily as Ninv = 1 [25] . Finally, and most importantly, we 
already know that SM neutrinos exist and behave like invisible particles at colliders. 
The decay chains in Fig. 1 can easily^ contain SM neutrino particles, which would 
contribute to the total invisible particle count Nin^ [26]. Under those circumstances, 
one needs to generalize the Mt2 variable to account for the extra invisible particles, 
and define corresponding variables Mt3, Mta, etc. [4], where the numerical subscript 
indicates the number of invisible particles Ninv In doing so, however, one faces the 
following fundamental problem. Recall that the most useful property of the original 
M'T2 variable was that its endpoint M^"^ equals the true parent mass Ma when the 
test child mass Mc coincides with the true child mass Mc- 

M^riMc,=McJ = MA- (1.2) 

However, as we shall explicitly see below in Section 3, this property depends crucially 
on the assumption that there aren't any additional invisible particles in the game. 

^However, see [23] for an update to [22]. 

^The standard example is the chargino decay in supersymmetry xt ~^ ^^^^Xii which gives two invisible 
particles — a neutrino u and a neutralino Xi- 
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In reality, whenever the cohection of particles V in Fig. 1(b) contains other invisible 
particles x addition to the invisible child particle C, and there are intermediate on- 
shell resonances in the cascade decay chain, the bound of eq. (1.2) is not necessarily 
saturated, and one can only write 



The inequality applies even when the true masses and Mc are being used. The 
same holds for the new "generalized" variables Mt3, Mt4, etc. from [4]. In fact, 
these generalized variables are simply related to the asymmetric M'r2 variable [19] as 
follows"^ 



and so on. Therefore, the endpoints of Mt3, Mta^ etc. are also not guaranteed to 
provide a saturated bound like eq. (1.2) and instead the best one can do with them 
is to put a lower limit on Ma as in eq. (1.3). One of the main goals of this paper, 
therefore, will be to address the problem of multiple invisible particles and propose 
how to recover saturated bounds of the type in eq. (1.2). 

(f) The final issue is inherently related to the previous point: once we allow for multiple 
invisible particles (A^i™ > 1) in the decay chains, we must also address the question 
of the correct event topology. In other words, we must ask the question, which 
Feynman diagram is hiding behind the yellow-shaded placeholder in Fig. 1(a). Most 
studies in the literature already assume that the correct event topology is known and 
rarely discuss what happens when this assumption is incorrect [27-31]. Consider, for 
example, the simplest possible case of one visible particle on each side {Nms = N^-^ = 
1) and then let us allow up to two invisible particles in a given decay chain {Ninv < 2 
and Nl^^ < 2). Even for this simple case, there are 10 possible decay topologies which 
are explicitly shown in Fig. 2(a-j). Furthermore, since one cannot be absolutely sure 
that the two visible particles originated from opposite decay chains, in principle one 
should also contemplate event topologies with = 2 and N^^^ = 0, and two such 
examples are shown in Fig. 2(k-l). 

The main goal of this paper is to illustrate the application of the Oxbridge M'r2 variable 
in the more general situations described in (c-f) above. We shall be particularly concerned 
with the issue (e) of multiple invisible particles and the related problem (f) of the unknown 
decay topology. It is quite difficult to construct a kinematic variable (or more generally, 
some kind of object) which would single out the correct event topology. An initial attempt 
in this direction was made in Ref. [32] which considered the decay of a single resonance A 
(as opposed to pair production as in Fig. 1) and studied the invariant mass distribution 
My of the visible decay products. It was suggested that there exists a correlation between 

^Here and below our notation is that the additional invisible particle Xi appears in the same decay chain 
as the invisible child particle d. See Fig. 2 for explicit examples. 



MP2^''{M^ + Mc) < Ma. 



(1.3) 



Mrs {Mc, , ■,Mc,) = Mt2 (M^ + M^, ,Mc,), 
MTi{Mc, , ; Mc^ , ^X2 ) = {Mc, + M^i - Mc^ + ^X2 ) > 



(1.4) 
(1.5) 
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Figure 2. Possible event topologies with two visible particles vi and V2 in the final state and up 
to 2 invisible particles {Ci^Xi} P^r decay chain. Bi and B2 are intermediate on-shell resonances: 
Ma > Mb, > Mc, ■ 



the peak, the curvature at the peak, and the endpoint of the My distribution, which in 
favorable scenarios can be used to infer the decay topology and the number of invisibles 
Ninv However, this method requires Nms > 2 and cannot be used to discriminate the 
topologies of Fig. 2(a-j), for which Nyis = 1 and My = My. = const. 

It appears, therefore, that we need to study event topologies on a case by case basis, 
i.e., assume in turn each one of the event topologies from Fig. 2, then test for consistency 
with the data. Note that the event topologies differ from each other in two aspects: 

1. The number of invisible particles Niny. For example, the diagrams in Fig. 2(a,k,l) 
have Ninv = N-^^ = 1, the diagrams in Fig. 2(e-j) have N^nv = N'-^^ = 2, while the 
diagrams in Fig. 2(b-d) have Nmy = 2 and N^^^ = 1. 

2. Different number of intermediate on-shell resonances Bi. For example, the diagrams 
in Fig. 2(a,d,j,l) have no intermediate resonances, the diagrams in Fig. 2(b,c,g,i,k) 
have one, while the diagrams in Fig. 2(e,f,h) have two. 

As already discussed earlier, both of these effects jeopardize the conventional Mt2 approach, 
and some modifications are required. In principle, there are two possible solutions. The 
first, presented in Section 2, is to modify the conventional Mt2 definition on a case by 
case basis — depending on the assumed event topology. The new definitions will use the 
appropriate number of invisible particles, and will also take into account the relevant on- 
shell kinematic constraints — one constraint for each intermediate resonance B^. Because 
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the mass shell constraints are 3 + 1-dimensional relations, this approach does not use the 
transverse Mt2 variable per se, but the alternative 3+1 dimensional invariant mass variable 
M2 from [2], supplemented with the appropriate mass-shell constraints (see also [33, 34] 
and the analogous constrained variables M2C and Msc). The advantage of this approach is 
that we do not have to change the interpretation of the corresponding kinematic endpoint 
jy^max — g|.jj| provides a saturated bound on the parent mass in terms of the masses of 
the invisible daughter particles, similarly to eq. (1.2). 

On the other hand, from a practical standpoint, the constrained M2 method has a 
certain disadvantage: one has to define a new M2 variable for each assumed event topology. 
In the majority of the cases, those variables are not equivalent to the conventional Mt2, 
and therefore cannot be calculated with the existing Mt2 codes [21, 22]. Since analytical 
formulas are also unavailable, one would have to write a separate code for each variable. 
The mass-shell constraints represent an additional complication, since most likely they 
will have to be solved analytically first [35-37]. This is why in Section 3 we present an 
alternative approach whereby one keeps the original definition of Mt2 (thus being able 
to recycle the existing numerical codes), and instead modifies the interpretation of the 
kinematic endpoint M^2^ on a case by case basis. The basic idea was raised in [38], where 
it was applied to the Mct variable [39] in the example of chargino decay. Here we show that 
one can apply a similar treatment to the Mt2 kinematic endpoints, and we also provide 
the recipe for a general event topology. 

In Section 4 and Section 5 we extend our method to the case of different children 
particles (Ci and C2) and diff^erent parent particles {Ai and A2), respectively. Finally, in 
Section 6 we study the impact of the number of invisible particles N^v on the shape of 
the Mt2^ distribution [40] . The technique was originally proposed in [26], which focused 
on the shape of the Mj'2 variable. Here we prefer to consider the doubly projected version 
Mj^2x) which is less sensitive to extraneous factors like spin correlations, the underlying 
partonic CM energy -y/s, etc. [41]. Thus our analysis completes the generalization of Mt2 
along the lines (c-f) discussed above. Since the directions (a) and (b) have already been 
discussed extensively in the literature, for the sake of clarity and simplicity, throughout 
the paper we shall use the simplest assumption N^ig = N'^^g = 1 and we shall introduce 
upstream Pt (i.e., transverse momentum for the parent system (A, A)) only as necessary. 

2 Maximally constrained invariant mass variables 

In this section we revisit the event topologies from Fig. 2 and for each one define the 
appropriate invariant mass variable which provides the maximal bound on the parent mass 
Ma [2]. 

The event topology of Fig. 2(a). This is the classic Mt2 event topology from Fig. 1(b), 
so the relevant variable is simply Mt2 [3] (or its asymmetric version if the children have 
different masses [19]). For future reference, it is convenient to introduce here the 3 + 1- 
dimensional version of Mt2, denoted simply as M2^, which will make it very easy to 
incorporate additional mass-shell constraints later on. 

^Not to be confused with the wino mass parameter in supersymmetry. 
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M2 is defined on a case by case basis, following the general recipe outlined in [2]. 
The particles in the final state are divided into two groups® (hence the subscript "2" on 
M2), one for each parent. Then, the larger of the two parent masses is minimized over 
the momenta of the invisible particles. The minimization is performed subject to all ex- 
isting kinematic constraints and assumptions. The measurement of the missing transverse 
momentum always provides one constraint on the net sum of the invisible transverse 
momenta, and another constraint is due to the assumption that the parents have equal 
masses (see discussion of item (d) in the Introduction) . 

In the case of the diagram in Fig. 2(a), there are no further constraints, and we get 

Ml^^-^iMc,;Mc,) = min {{P,^ + Pc,f} (2.1) 



{P,^+Pc^f = {P,^+Pc, 



Ptc^ + PtC2 = fT 

Upon performing the minimization over the longitudinal momenta PzC\ find PzC^^ one finds 
that eq. (2.1) is equivalent [2, 33] to the usual Mt2 when the children's masses are taken to 
be the same 

M2ia){Mc; Mc) = Mt2{Mc), (2.2) 

or to the asymmetric M'T2 [19] 

M2(,) (Mc, ; Mc72 ) = Mt2 {Mc, ^Mc^), (2.3) 

if the children's masses are kept different. 

The event topology of Fig. 2(b). Proceeding as before, for the diagram in Fig. 2(b) we 

get 

Ml^^{M^,,Mc,:Mc,)= min {{P^^ + P^^ + P^^f] . (2.4) 

^Xl '■'iCi ,^C2 



{P,^+P^^+Pc,f = iP,,+Pa, 



Pi = Ml 

Pc, = 

Pk = 
{Px.+Pc,? = Ml, 

Ptxi + Ptc, + Ptc2 = fT 
With Pb, = -Pxi + Pc\ ) this can be equivalently written as 

M2%) {M^, , Mc, ■,Mc,) = ^min { {P,, + Pb,?] ■ (2.5) 



PB^,Pa 



2 



®The division into two groups is motivated by the two parent hypothesis. If we do not make this 
assumption and instead treat the event as a whole, we are led to the global inclusive variable VSmin [42, 43], 
which in the language of [2] is denoted as M\. 
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Pk = 

Pc, = M'c, 
Ptb, + Ptc2 = 

From here, comparing to eq. (2.1) and taking into account eq. (2.3), it follows that 

M2(b) (M^, , Mc\ ]Mc,) = M2(,) {Mb, ; Mc^ ) = Mt2 {Mb, ,Mc,). (2.6) 

As expected, the kinematics of Fig. 2(b) is described by the asymmetric Mt2 variable [19], 
with the intermediate invisible resonance Bi treated effectively as a final state invisible 
particle. 

The event topology of Fig. 2(c). We get 

mI(^^){M^„Mc,;Mc,) = ^ min {{P^^ + P^^ + Pc^f} . (2.7) 



{P,^+P^^+Pc,f = {P., + Pc, 



\2 



Pi = Ml 
Pc, = Mk 
Pk = Mk 
{P,,+Pc,f = Mk 
Ptxi + Ptc, + Ptc2 = fr 

This variable is new — it cannot be related to existing versions of the Mt2 variables as in 
eqs. (2.3, 2.6). 

The event topology of Fig. 2(d). Here there are no mass-shell constraints and we have 
^l(d)(^xi'^C,;Mc,) = min {(P,^+P^^+P^J2|_ (2.8) 

(P,,+P^,+PCi)' = {PV2^PC2? 

pk = Mk 
pk = Mk 

Pk = Mk 
Ptxi + Ptc, + Ptc2 = fr 

It is easy to show that the minimization in eq. (2.8) selects the collinear momentum con- 
figuration [2, 42] 

P =M2ap (2.9) 

and eq. (2.8) can be rewritten as 

M\^.^{M^„Mc,:Mc,) = min {{P,^+P^^f] . (2.10) 

(P., +P*J2 ^ (p^^+p^j2 
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Ptv&i + Ptc2 = fr 
in terms of an effective composite invisible particle ^'i with mass 

M^,=M^,+Mc, (2.11) 

and 3-momentum 

P,,^P,,+Pc, = ^Pc. = ^P,,. (2.12) 
Comparing to eq. (2.1), we recognize this as the asymmetric Mt2 variable in eq. (2.3) [19] 
M2(rf) (M^, , Mc, ■,Mc.,) = M2^a) {My,, + Mc, ■,Mc,) = Mt2 (M^i + ^c, .Mo,). (2.13) 
The event topology of Fig. 2(e). This is our first example with four invisible particles: 
Ml,^{M^,,Mc,;M^,,Mc,)= min [[p^^ + + p^^f] . (2.14) 

^Xl '-M^i :-rx2 '^C2 

(P., + P^, + Pc,? = (P.2 + Px2 + Pc,? 

Pi = 

Pc, = Ml^ 
Pi = 

Ph = 
{P^,+Pc,? = Ml^ 

{P^, + Pc,? = Ml^ 
Ptxi + Prcr + Ptx2 + Ptc2 = fr 

By introducing Pg. = P^. + Pc-, {i = 1,2), this can be equivalently rewritten as 

M\^^{M^,,Mc,-,M^,,Mc,) = ^uvin {{P,^+PB^f] (2.15) 

\2 



{P,,+Pb,Y = {P,2+Pb2Y 

Pk = Ml 

Pk = Ml 
Ptb, + Ptb2 = fr 

and therefore is again reduced to the asymmetric Mt2 [19] 

M2(e) (Mxi > Mc, ; Mx2 ,Mc2) = M2(a) {Mb, ■,Mb2) = Mt2 {Mb, .Mb,)- (2.16) 

Notice that M2(e) does not depend on the test masses M^i' Mc,, M^, and Mc,- In fact, 
the two invisible decays Bi ^ Xi + Ci, (i = 1, 2), have no observable consequences. 
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The event topology of Fig. 2(f). This is similar to the previous case, except the order 
in which particles xi ^-^id vi appear in the decay chain is reversed 

Mlt^f)iMxi^Mc,-,Mx2,Mc,) = ^ min {(P^^ + P^^ + P^^f] . (2.17) 

^Xl 1-^X2 '^^2 



(P., + P^, + Pc,? = {Pv, + Px2 + Pc,f 



Pi = Ml 

Pc, = 
Pi = Ml 

Pk = 

{P.,+Pc,f = Ml^ 

{Px. + Pc,? = Ml, 
Ptxi + Ptc^ + Ptx2 + Ptc, = fT 

Once again, the B2 mass shell constraint causes the inputs My^^ and Mc^ to drop out and 
we get 

Ml^^^iM^„Mc^;Mx2,Mc,) = ^ min {{P,^ + P^^ + P^^f} . (2.18) 

-fxi >-^Ci i-f^2 



(P., = {Pv,+Pb.. 



2)' 



Pi = Ml 
Pk = Ml, 
PI2 = Ml, 
{P,,+Pc,f = Ml, 

Ptxi + Ptc^ + Ptb2 = 
This variable is also different from Mt2, but is related to M2(c): 

M2(/) {M^, , Mc, ; Mx2 ,Mc,) = M2(e) (M^i ' Mc, ■,Mb,). (2.19) 

The event topology of Fig. 2(g). This case is similar to Fig. 2(e), but the Bi mass shell 
constraint is removed: 

M^^g)iMxi,Mc,;Mx2,Mc,) = ^ min {(P^^ + P^^ + P^^f} . (2.20) 

^Xl >-fX'i '-fx2 '^C2 
(P„, + P^, + Pc,f = (P.2 + PX2 + ^^02)' 

^xi = Ml 
Pc, = Ml, 
Pi = Ml 

pk = Ml, 

2 11^-2 



{Px2+Pc,r = Ms^ 

Ptxi + Ptc^ + Ptx2 + Pre, = fr 
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Introducing the effective particle ^ as in eqs. (2.11, 2.12), we get 

M2%)(M^„Mci;Mx2,Mc,) = min {{P,^+P^^f} . (2.21) 

2 f tCj < T\~J \2 



P4^ = (Mx.+Mc, 
Pk = M 
Pt^, + PtB2 = fT 



and tfierefore 

M2(g) (M^, , Mci ; > ) = M2(„) (M^, + Mc, ; Mb J = Mt2 {Mxi + , M^, ) . (2.22) 
The event topology of Fig. 2(h). This is another new and non-trivial case: 

Ml^h)iMxi,Mc,;Mx„Mc,) = ^ min {(P^^ + P^^ + P^^f} . (2.23) 



(P., + P^, + Pc,? = {Pv, + Px2 + ^'C: 



■J 



Pi = Ml 

Pk = Ml, 

PL = Ml 
Ph = 

{P.,+Pc,? = Ml^ 
{P,, + Pc,? = Ml^ 
Ptxi + Ptc^ + Ptx2 + Pre, = 

The variable M2(^h) cannot be reduced to one of the previous variables and would have to 
be evaluated separately. 

An interesting variation of the M2{^h) variable arises in the symmetric case when the 
intermediate on-shell particles Bi and B2 are the same, with = M^j. Then, one can 
replace the two mass shell constraints for Bi and B2 with the requirement that the Bi and 
B2 masses are equal, but without specifying the actual numerical value: 

M\h)iMx,,Mc,;Mx2^Mc,) = ^ min {[p^^ + p^^ + p^^f] . (2.24) 

^Xl >-ft7i '-fx2 '^C2 
(P„, + P^, + Pc,f = (P.2 + PX2 + PC2? 

Pi = Ml 
Pc, = Ml;, 
Pi = Ml 

Pk = Mh, 

{Pv,+Pc,? = {Pv,+Pc,? 
Ptxi + Ptc^ + ^¥x2 + Ptc2 = h 
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The advantage of this approach is that one does not need to know the value of Mb^ = 
beforehand. 

The event topology of Fig. 2(i). This is similar to the previous case, except the B2 mass 
shell constraint is absent: 

M|(,)(M^„Mc,;Mx2,^c.) = ^ min {(p^^ + + p^.^f} . (2.25) 



(Pv, + Px. + Pc.? = {Pv, + Px, + PC: 



Px. = ^l. 

Pk = 

Px. = ^X. 

Pc, = 



Ptxi + Ptc^ + Ptx2 + Ptc2 = fT 
With the help of the effective invisible particle this can be again simplified: 

M%){Mxi^Mc,;M^,,Mc,)= min { (P,^ + + _ (2.26) 

2 _ / D I D \2 



p,^ + p^^ + Pc,Y = {p,, + p^,) 



Pi = Ml 
Pc, = 

Pi, = {M^,+Mc,? 
{P.,+Pc,? = Ml^ 

Ptxi + PtCi + Pt^2 = fT 
Comparing to eq. (2.7), we see that this variable is related to M2(c)- 

M2(i) (M^, , Mci ; Mx2 ,Mc2) = M2(e) [M^i > ; + Mc,). (2.27) 
The event topology of Fig. 2(j). Here there are no mass shell constraints and we have 
M|(^.)(M^„Mc,;M^2,Mc,) = min {(P^^ + P^^ + P^^f} . (2.28) 

2 _ / D , D , D„ \2 



(P., + p^, + Pc,r = (P.2 + Px2 + PC2) 

Pi = Ml 
Pk = 
Pi = Ml 
Pc2 = Ml, 
Ptxi + Ptc^ + Ptx2 + Ptc2 = h 

Introducing two effective invisible particles ^'i and 'J/2, this reduces to 

M^.^^{M^„McV,Mx2,Mc,) = min {(P,^ +P^j2| _ (2.29) 
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Pi = (Mx.+Mcf 
Pi, = {M^,+Mc,? 
-Pr^i + Pt^2 = 

Comparing to eq. (2.1), we see that 

M2(j) (M^i , Mci ; ,Mc,) = M2(a) (Mxi +^c^ ; M^, +Mc,) = Mt2 [M^.+Mc^, +^0, 

(2.30) 

In conclusion of this section, let us summarize its main points. We considered the 
ten event topologies in Fig. 2(a-j) and defined the corresponding maximally constrained 
invariant mass variables, which fell into three categories. 

• The variable M2(a) and its cousins M2{h)i -A^2(d); -^2(e); -^2(g)i fiiid M2(j), all of which 
can be computed in terms of the asymmetric Mt2- 

M2^b) {Mx, , Mc, ; McJ = M2(„) {Mb, ; Mc, ) 

= Mt2{Mb„Mc2) 

M2(d) (Mxi > ^c, ■,Mc2) = ^2(a) {My,, + Mc, ; Mc, ) 

= Mt2{Mxi+Mc„Mc2) 

M2(e) {Mxi ' Mc, ; Mx2 ,Mc, ) = M2(„) {Mb, ; Mb, ) 

= Mt2{Mb„Mb,) 

M2(g) {M^, , Mc\ ; Mx2 ,Mc,) = M2(„) {M^, + Mc, ; Mb, ) 

= Mt2{Mxi+Mc„Mb,) 

M2(,) {M^, , Mc\ ; Mx2 , ) = ^2(a) (M^i + ; Mx2 + Mc, ) 

= Mt2(Mxi + Mc, , Mx2 + Mc,) 

• The variable M2(c) and its friends M2(^f) and M2(j): 

M2(y) (M^i , Mci ; M^2 ,Mc,) = M2(e) (Mxi > ; Mb, ) 
M2(,) (M^i > ; Mx2 ,Mc,) = M2(e) (Mxi > ^c, ; M^, + Mc, ) 

• The variable M2{h)- 

The advantage of these maximally constrained invariant mass variables is that their end- 
points provide saturated bounds on the parent mass Ma- In other words, we recover the 
main feature in eq. (1.2) of Mt2 which was lost due to the presence of the mass shell con- 
straints. More specifically, when we choose the test masses M^. and to be the true 
ones (denoted without a tilde) , the corresponding M2 kinematic endpoint gives the parent 
mass Ma'- 

M^^^^^{M^^,McJ = Ma. (2.31) 
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Mt2 calculation 
V 



A 



C 



A 



C 



V 

(a) 



MJpf^ interpretation 



A 



V 



^1 



A 



*5 



V 



(b) 



Figure 3. The effective topology used for (a) the calculation of the Mt2 variable and (b) the 
interpretation of its endpoint M™2°^ in the different cases from Fig. 2. In panel (b), ^"1 and ^'2 are 
effective particles whose masses are calculated in terms of the physical masses M^, M^, Mc and 
according to the rules in Fig. 4. 

Here the subscript (...) stands for the event topology specifier: (a), (6), ....(j). This also 
reveals the main problem with the maximally constrained variables M2(..,) — while we re- 
covered a saturated bound in eq. (2.31), we only managed to do so at the cost of introducing 
a separate M2 variable for each event topology. Furthermore, in many cases, these variables 
cannot be calculated with the existing public codes. This motivates an alternative, more 
practical approach, which will be the subject of the next section. 

3 Effective event topology and reinterpretation of the kinematic end- 
point of the usual Mt2 variable 

3.1 General setup 

In this section, we revisit the 10 event topologies in Fig. 2(a-j) and insist that we analyze 
all of them by means of the conventional Mt2 variable, which is computable by the public 
codes. In other words, for the calculation of the Mj'2 variable, we will disregard any 
differences between the diagrams in Fig. 2(a-j) and instead pretend that the events arise 
from the classic Mt2 event topology shown in Fig. 3(a): two parents with equal masses Ma 
are produced, and each one subsequently undergoes a two body decay to a single visible 
particle V with constant mass'' My = const and a single invisible particle C with mass 
Mc- As usual, we assume that the two children particles are identical (or at the very least, 
that they have a common mass Mc)- 

Following the usual procedure, we can then build the Mt2{Mc) distribution and ex- 
tract its endpoint M^"^(Mc). It is only at this point that we need to worry about the 
actual origin of the events. In order to account for the differences between the diagrams 
in Fig. 2(a-j), we propose to interpret the measured endpoint M^^^ {Mc) in terms of the 

^For simplicity, in this paper we consider massless V particles, i.e. My ~ 0. 
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A 



C 



= Mc 
, X V 



B 



C 



A/. = Af^{l-A^(l-||)e"} 



1/2 



id) 



B I . 

Af^ = Mb 
v\ ,x 



C 



C 



A'U = M^+Mc 



Figure 4. Equivalence diagrams for the interpretation of M^°^ in Fig 3(b). Each decay chain 
(a-d) can be replaced with the effective vertex (e). The mass of the effective particle 5* should 
be suitably chosen as shown in each panel (a-d). 



effective event topology shown in Fig. 3(b): we still produce two identical parents, but now 
they decay to effective invisible particles ^'i and ^'2, respectively. The interpretation of 
the endpoint M^2^{Mc) is still given by the usual formula [15, 16] 



{Mi 



c) 



(3.1) 



where 



Ma 
2 



Ml 



Ml 



1/2 



(3.2) 



and it is only the masses M^^ and Mq;^ of the effective particles that would have to be 
calculated on a case by case basis, depending on the assumed event topology from Fig. 2. 
The effective masses M^-,{i = 1,2) are in general functions of the true masses of the 
particles appearing in the corresponding event topology from Fig. 2. In order to derive those 
functions, we use the technique originally suggested in [38] for the case of the contransverse 
variable Mqt- The main idea is to consider the extreme momentum configuration which 
gives the maximum value M^"^. The rules for constructing the effective mass M,j, for a 
given cascade decay chain are listed in Figure 4. Most of them should be pretty intuitive, 
given our discussion in the previous section, where we already had to introduce an effective 
particle ^. For example, we already encountered Fig. 4(d) in eq. (2.11) and we also saw 
the replacement of Fig. 4(b) in eq. (2.6). Perhaps the one non-trivial example is that of 
Fig. 4(c), where the effective particle mass is formed as 



M^, {Ma, Mb, , Mc, , M^J = Ma{i 



Mb, 
Ma 



Ml 



1/2 



(3.3) 



where 



cosh"^ 



m1-fM 



Ml 



2MaMb, 



(3.4) 
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When the invisible particle x is massless (M^ 
reduces to 

M^, (Ma, Mb, , Mc, , 0) = Mc, 



0), the effective mass M^- in eq. (3.3) 

Ma 



Mb, 



(3.5) 



The advantage of the Mt2 endpoint reinterpretation approach is that one never has 
to stray from the conventional procedure — it is only at the very last stage of the analysis 
that the event topology issue comes into play. The M^2^ endpoint is always interpreted 
universally as in eq. (3.1), and only the parameter n from eq. (3.2) depends on the event 
topology from Fig. 2: 



M(e) 



Ha) 



Ma 
2 

Ma 
2 

Ma 
2 

Ma 
2 



1 



Ml 



1 



Mh 
Ml 



\ 




Ml 



Ml 



Ma 



Mb, 
\\ Ma 



Ml, 




Ml, 



Mh 
Ml 



\ 



{Mc,+M^,Y' 



Ml 



Ml 



\ 




M' 



B2 



Ml 



Ml 



Ma 



Mb, 
^ Ma 



Ml 



Ml 



Ma 
2 



\ 



{Mc,+M^,] 



Ml 



Ml 



Ma 



Mb, Mb, 
Ml 



Ml 



M. 



Ml 



em/2 
Ma^ 



Ma 



Mb, 
\\ Ma 



Ml 



{Mc, + M- 



X2J 



Ml 



1 



{Mc,+M^,] 
Ml 



1 



{Mc,+M^,f 
Ml 



(3.6) 

(3.7) 

(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 

(3.14) 
(3.15) 



with r}i still given by eq. (3.4). 

3.2 Application to chargino decays 

A well-motivated and relevant physics case illustrating the decay topologies in Fig. 4 is 
provided by the chargino decays in supersymmetry. If we identify the parent particle A 
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with the chargino and the visible SM particle with a charged lepton the invisible 
daughter particles could be: a neutralino x^, a sneutrino z>, or a SM neutrino u. Then each 
of the decay chains in Fig. 4(a-d) can be interpreted as follows: 

• The topology of Fig. 4(a): A ^ v + C. This corresponds to a two-body decay of 
the chargino to a lepton and sneutrino: i'^ + with the sneutrino being the 
lightest supersymmetric particle (LSP). 

• The topology of Fig. 4(b): A ^ v + B, followed by S — )■ x + C- This is similar to the 
previous case, where we identify B with the sneutrino: x~^ + D, only this time 
the sneutrino decays further invisibly: D ^ u + x^. 

• The topology of Fig. 4(c): A — )■ x + followed hy B ^ v + C. Here there are 
two possible examples. The first one is analogous to the sneutrino decay considered 
above — only this time the intermediate particle is a charged slepton and we get 
X~^ — )■ v + £~^, followed by — t- + x^ • The invisible particle masses are identified 
as Mc = My.0 and = My = 0. Another possibility is to have the chargino decay 
to an on-shell W boson: x^ + , followed by + v, in which case 
the invisible masses are Mq = My = and M^ = M^o. 

• The topology of Fig. 4(d): ^ — t- f + x + C- This is realized if the chargino two-body 
decays are closed and it decays via a virtual slepton or W boson: x"*" + 1^ + x^- 

By pairing up these 4 cases, we can obtain all 10 event topologies from Fig. 2(a-j). 

As a specific scenario realizing these patterns, let us consider the Tchislepslep sim- 
plified model [44] , where > Mj>^ ~ M^+ > M^o . The masses of the sneutrino and 
the charged slepton are taken to be equal because they belong to the same SU{2)w dou- 
blet (they are both left-handed). The two possible chargino decays are given by Fig. 4(b) 
and Fig. 4(c) and have equal branching fractions. The three possible ways of pairing up 
Fig. 4(b) and Fig. 4(c) lead to the event topologies of Fig. 2(e), Fig. 2(f) and Fig. 2(h): 

• Sneutrino- sneutrino. When both charginos decay through a sneutrino as in Fig. 4(b), 
we obtain the event topology of Fig. 2(e) with the identifications Ma = -^x+' -^^i ~ 

= My, Mc^ = Mc^ = M^o, and M^^ = M^^ = My = 0. Therefore the 
Mt2 endpoint in eq. (3.1) should be interpreted with the /i(e) parameter taken from 
eq. (3.10): 

f M?\ , , 

m = /^^^ = ^ 1 1 - I • (3-16) 



Slepton- sneutrino. The hybrid pairing of Fig. 4(b) and Fig. 4(c) leads to the event 
topology of Fig. 2(f) with the identifications Ma = M^+, Mb^ = M^+, Mb^ = My, 
Mci = Mc2 = ^x"' ^xi ~ — — 0- The Mt2 endpoint for such events 
is interpreted with ^(j) from eq. (3.11): 



\ 




(3.17) 
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Figure 5. Distributions of the Mt2 variable (for zero test mass Mc = 0) for the Tchislepslep 
simphfied model, with = 500 GeV, M^,^ ~ M^-+ = 400 GeV and M^o = 100 GeV. The chargino 

branching fractions are fixed as B{x^ l'^ + v) = B{x^ v + — 50%. Results are shown 
for direct chargino pair production at the LHC with 8 TeV CM energy. The vertical dashed lines 
mark the expected endpoints for each event topology: sneutrino-sneutrino (red), sneutrino-slepton 
(blue) and slepton-slepton (green). 



• Slepton-slepton. When both charginos decay through a slepton as in Fig. 4(c), we 
obtain the event topology of Fig. 2(h) with the identifications Ma = Mb^ = 

= M^+, Mc^ = = M^o, and M^^ = M^^ = M^ = 0. The Mt2 endpoint for 
such events is interpreted with from eq. (3.13): 




(3.18) 



The Tchislepslep scenario is illustrated in Figure 5, where we show Mt2 distributions 
from Monte Carlo simulations for direct chargino pair production, at the LHC with 8 
TeV CM energy. The electroweak mass spectrum is chosen as = 500 GeV, Mc,j^ — 
Mj;+ = 400 GeV and M^o = 100 GeV, while the chargino branching fractions are fixed as 

B{x~^ — 7- ^+ + P) = B(x~^ V + i'^) = 50%. For simplicity, we calculate Mt2 with zero 
test mass {Mc = 0), so that the endpoint formula in eq. (3.1) simplifies to 

180 GeV for sneutrino — sneutrino events, 
290 GeV for slepton — sneutrino events, (3.19) 
469 GeV for slepton — slepton events. 



M^^^'iMc = 0) = 2// = < 
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where the /i parameter was calculated from eqs. (3.16-3.18), depending on the event type. 
In Figure 5, sneutrino-sneutrino events are plotted in red, slepton-sneutrino events are 
plotted in blue and slepton-slepton events are plotted in green. The color-coded vertical 
dashed lines mark the corresponding Mt2 endpoints expected from eqs. (3.16-3.19). 

We conclude this section with a discussion of the TChiww SMS model, in which the 
sleptons and sneutrinos are heavy, and both charginos decay through an on-shell W boson 
instead (provided, of course, that > M\y + M^o). The corresponding topology is the 
one from Fig. 2(h), where Ci = v and Xi = x'^ is the neutralino LSP. The effective masses 
from eq. (3.3) are 



where 



rj = cosh" 




Mw 




2 

W 



Mir 



2M^+Mw 



Then the Mt2 endpoint (with zero test mass M^o) is given by 

{0) = Mwe'^. 



(3.20) 



(3.21) 



(3.22) 



Finally, if the mass splitting between the chargino and the neutralino is small iM^+ — 
M^o < Mw), the chargino decays as in Fig. 4(d). The event topology is given in Fig. 2(j) 
and the Mt2 endpoint (again with zero test mass M^o) is 



M^s""^ (0) 



Ago 



(3.23) 



4 Reinterpretation of the kinematic endpoint of the asymmetric Mt2 
variable 



In the effective topology of Fig. 3(a), the test masses for the two invisible children were 
chosen to be the same: Mq- However, as demonstrated in [19], it is straightforward to 
generalize the standard Mt2 calculation in Fig. 3(a) to allow for different test masses for 
the two children particles. In other words, we could use the effective diagram in Fig. 3(b) 
not only for the interpretation of the kinematic endpoint, but also for the actual calculation 
of the (asymmetric) Mt2 variable itself, in terms of two different test masses'^, M^^^ and 
M^^ [19]: 



mill {max{MTAi(^'T*i,Af*i), Mta2(^t*2.^^*2)}} : 



(4.1) 



where MxAi is the transverse mass of the parent particle A- 

MtaAPt9,,M9,) 



Ml +Ml^+2 



(4.2) 



^Recall our notation that test masses always carry a tilde. 
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and the subscript "D" in eq. (4.1) is used to remind the reader that the asymmetric Mt2 
variable uses two different test mass inputs. In eq. (4.2), My. and PrVi are correspondingly 
the invariant mass and transverse momentum of the effective visible particle Vi resulting 
from the decay of Ai. Similarly, Mij,. and -Pripj are the test mass and test transverse 
momentum of the effective invisible particle ^'j (see Fig. 3(b)). As usual, -Pt*i and -Pt*2 
are subject to the missing transverse momentum constraint, and then varied to find the 
minimum of the function in eq. (4.1). 

When the visible particles are massless (My- = 0), as we are considering here, the 
kinematic endpoint of Mt2D is given by [19] 



M^2Z)(^^*i-^*2) 




1/2 



(4.3) 



where the parameter n encoding the dependence on the physical masses is still given by 
eq. (3.2), while the two test masses enter through the combinations 

Ml = ^(m|^+M|J, (4.4) 



M 



2 



1 



M|^ - M|^ 



(4.5) 



Obviously, one can always go from the asymmetric Mt2D variable in eq. (4.1) back to 
the original Mt2 variable [3], simply by setting the two test masses to be equal: 

Mt2d{M^M^) = MT2(Af*). (4.6) 

Furthermore, the Mt2D endpoint in eq. (4.3) only allows us to measure one parameter: //, 
and the same can be said about the usual Mt2 endpoint in eq. (3.1) as well. A natural 
question then is whether there is any benefit at all from introducing the more complicated 
variable Mt2D- We see two motivations for considering Mt2D- 

• In the context of manifestly asymmetric event topologies like the ones shown in 
Fig. 2(b), Fig. 2(c), Fig. 2(d), Fig. 2(f), Fig. 2(g) and Fig. 2(i), the language of Mt2D 
is more appropriate because the masses of the effective invisible particles ^'i and ^2 
are different, even if the masses of the particles Ci and C2 at the end of the decay 
chains are the same. 

• More importantly, by considering the asymmetric Mt2D variable, one could in princi- 
ple find not just one constraint among the three unknown parameters M^, Mijr^ and 
Mijij, but the actual values of all three parameters themselves [19, 20]. To this end, 
one needs to consider events in which the parent AA system recoils against upstream 
objects (from initial state radiation or prior decays) with a net upstream transverse 
momentum Ut- One then measures the endpoint in eq. (4.3) and compares the results 
in different Ut bins. In general, the results for the Mt2D endpoint will depend on 
the value of Ut- However, when the test masses Mij,. are equal to the corresponding 
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Figure 6. (a) A generic event topology with different parent particles and (b) the effective event 
topology for its interpretation as discussed in the text. 



true values Mij,^, the Ut dependence disappears and M^2d becomes independent of 
the upstream momentum [19] : 



0. 



(4.7) 



=M, 



The true location of the invisible masses M^- is often also revealed as the crossing 
point of several creases in the two-dimensional hyper-surface defined by the function 



[20]. 



5 Interpretation in the case of different parent particles 

As another application of the effective topology method, in this section we revisit the case 
of different parent particles {Ma-, ^^Ma^ as illustrated in Fig. 6(a). A well- motivated 
example of such a "coproduction" channel is provided by associated gluino-squark produc- 
tion which was studied in [24, 45], finding a correlation between the endpoint MJ^2^ of the 
conventional Mt2 variable and the larger of the two parent masses: 



{max) 



T2 



(Mc = Mc) = max {Ma, , Ma^ 



(5.1) 



This result, also suggested in [46], is the analogue of eq. (1.2) for the case of different 
parents. Here we point out that the relation in eq. (5.1) does not hold in general. The 
main result in this section will be the proper interpretation of the Mt2 kinematic endpoint 
in the case of different parents (Fig. 6(a)). The discussion will be organized as follows: in 
Section 5.1 we first treat the case with no upstream momentum {Ut = 0) and then in 
Section 5.2 we discuss the more general case with Ut 7^ 0. 

5.1 Events w^ith no upstream momentum {Ut = 0) 

We start from the known fact that the Mt2 endpoint is invariant under "back-to-back" 
boosts in the transverse plane [16]. In the conventional case of Fig. 1(b) with identical 
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parents, this boost brings both parent particles to rest. However, when the parents have 
different masses as in Fig. 6(a), no back-to-back boost can bring the parents to rest simulta- 
neously, and the "back-to-back" boost invariance of MJpf^ will be lost. More specifically, if 
the parent system (Ai, A2) is produced with some total CM energy v^, the parent particles 
Ai are boosted by the respective factors 



r?i(v^) = cosh 
^?2(^/^) = cosh 



1 / g+Mi^-Mi; 

\ 2s Ma, 

1 / g + Mi^-Mi; 

I 2s Ma^ 



(5.2) 
(5.3) 



For definiteness and without loss of generality, from now on we shall assume that Ma2 > 
Mai, so that max [Ma,, MA2) = -^^2- Then from eqs. (5.2, 5.3) it also follows that r]2 < 
Thus if we perform a "back-to-back" boost on the particles Ai and A2 with the boost factor 
r/2 corresponding to the heavier particle, the result will be that particle A2 will be brought 
to rest, while the lighter parent particle Ai will have some residual boost 



6r][Vs) = r/i(v^) - ri2{Vs) = cosh" 



Ml + Ml 



Ml - 



2Ma,Ma, 



(5.4) 



Comparing this to eq. (3.4), we see that eq. (3.4) and eq. (5.4) become identical if we identify 



Ml - Ml 



(5.5) 



The physical meaning of eq. (5.5) is the following — the events in Fig. 6(a), which represent 
pair production of different parent particles, are instead reinterpreted as in Fig. 6(b), which 
shows the pair production of identical parents A2, one of which decays into Ai plus a hypo- 
thetical invisible particle x whose mass is given by eq. (5.5). At this point it is important 
to note that the mass as defined in eq. (5.5) is not constant, but carries dependence on 
the CM energy y^, which is inherited from eqs. (5.2, 5.3). At hadron colliders, events are 
produced with varying thus in general takes values in the interval 



< M^(VI) < Ma, - Ma, , 



(5.6) 



where the lower bound corresponds to ^/S — )■ 00 and the upper bound is obtained at 
threshold VS = Ma, + Ma, ■ 

Once we have an event topology with equal parents as in Fig. 6(b), we already know 
how to interpret the Mt2 endpoint — we just follow the prescription from Section 3. The 
topology of Fig. 6(b) is of the type shown in Fig. 2(c), so for the upper decay chain we need 
to introduce an effective invisible particle ^'i as prescribed in Fig. 4(c) and eqs. (3.3, 3.4): 



Mvt, = MaAi 



Ma, 
Ma, 



1 



(5.7) 
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Then, the kinematic endpoint of the usual Mt2 variable is found from eq. (3.1): 



M^i^'iVs, Mc) = K^s) + V/u(VI)2 + M2 , (5.9) 
with /i(v^) given by eq. (3.8) 



\ Ma 




2 \ / A#2 



(5.10) 



Notice that the kinematic endpoint M^2^ in eq. (5.9) this time depends not only on 
the test child mass Mq, but also on the partonic CM energy At the LHC, ^/S is not 
constant, but varies from one event to another in accordance with the parton distribu- 
tion functions (PDFs). Therefore, the kinematic endpoint in eq. (5.9) will in general^ get 
smeared. As we can see from eqs. (5.4, 5.10), //(VI) is an increasing function of Since 

itself varies from its threshold value V§ = Mai + Ma2 to y/§ oo, the function /u(v^) 
takes values in 

fJ-min < //(v^) < fJ'max, (5-11) 



where 



flmin = hm n{Vs) 



Vs^Ma^+Ma^ 2 Ma. 

Ma, 




Ml 



(5.12) 



IJ-max = lim At(v^) 



Ml \ Ml\ IMa. 



J-^. (5.13) 



^ Mljy MlJ-^"''''\jMA, 



Now it follows from eq. (5.9) that even if we choose the test mass Mc to be the true mass 
Mc, the corresponding Mt2 kinematic endpoint M^pf^iyMc) will still vary with \/I between 
a minimum value of 



M^2""(Vs = Ma, + Ma,) = f^min + f^Ln + MS (5.14) 
and a maximum value of 

M^ri^S ^ oo) = f,ma. + ^Jfl'La.+Ml. (5.15) 

This is illustrated in Fig. 7, where we plot the Mx2 distributions for these two extreme 
cases: \/I = Ma-, + Ma, (blue) and — )• oo (red). We consider asymmetric events with 
different parents as in Fig. 6(a) and choose the mass spectrum as follows Ma, = 1000 GeV, 
Ma, = 200 GeV and Mc = 100 GeV. The blue and red vertical dashed lines mark the 



^Modulo the special case where Ai and A, are produced in the decay of some heavy narrow resonance 
X &s X ^ AiA, [41, 48]. 
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Figure 7. Unit-normalized Mt2 distributions for the asymmetric event topology of Fig. 6(a) with 
different parent particles Ai and A2 and no upstream momentum [Ut — 0). The particle mass 
spectrum is chosen as Ma^ = 1000 GeV, Ma^ = 200 GeV and Mc = 100 GeV. The Mt2 variable 
is computed with the correct value for the test mass Mc = Mc = 100 GeV. The blue histogram 
uses only events at threshold (a/I = Ma^ +MA2) and its expected endpoint, marked by the vertical 
blue dashed line, is given by eq. (5.16). The red histogram shows the corresponding result in the 
infinite energy limit \/i — > 00 (in practice, we take ^/~^ = lOOMyi^) and the expected endpoint 
from eq. (5.17) is denoted by the vertical red dashed line. The black vertical dashed line is the 
prediction from eq. (5.1). The green dotted line is the result from Pythia 6.4 [47] simulation at 
LHC8 for a realistic physics example corresponding to the event topology of Fig. 6(a): associated 
squark-chargino (q, x^) production, followed by g q + x° a^id Xi + v, for the same mass 

spectrum, (Mg, M^± , M^, M^o) = (1000,200,100,100) GeV. 

locations of the expected endpoints in eqs. (5.14) and (5.15), respectively. With the mass 
spectrum chosen for the figure, one gets 

M^^^'^i^ = Ma, + Ma2,Mc = Mc) = 409.8 GeV, (5.16) 
M^^^^'iV^ ^ 00, Mc = Mc) = 873.1 GeV. (5.17) 

Fig. 7 demonstrates that in these two limiting cases, the endpoints of the Mt2 distributions 
agree perfectly with our expectations in eqs. (5.16, 5.17) and stay well below the conjecture 
of eq. (5.1), which is indicated by the black vertical dashed line. 

For intermediate, more realistic values of VS, the upper endpoints of the corresponding 
Mt2 distributions will populate the region between those two extreme values, but will cer- 
tainly not exceed the theoretical maximum in eq. (5.17). As an illustration, in Fig. 7 we also 
show results (the green dotted histogram) from a realistic physics example simulated with 
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Pythia 6.4 [47]. We considered associated squark-chargino production in supersymmetry, 
pp — )• qx^, followed by g — >• q+x^ and xf ~^ i^+i^, with the same mass spectrum as before, 
(Mg, M^±, Mp, M^o) = (1000, 200, 100, 100) GeV. Such events fall into the different parent 
category of Fig. 6(a). With the available statistics, the Mt2 endpoint for the green dotted 
histogram happens to be around 800 GeV, which, as expected, is in between eq. (5.16) and 
eq. (5.17). We see that the realistic Mt2 distribution indeed does not saturate the bound of 
eq. (5.1). Therefore, the correct interpretation of the Mt2 endpoint in the case of different 
parents should be made with the help of eqs. (5.9, 5.10) instead. 

5.2 Events with upstream momentum (Ut 0) 

In the previous subsection 5.1, we considered events with no upstream momentum (JJt = 
0), where the parents Ai and A2 are produced back-to-back in the transverse plane. In 
reality, however, the inclusive (^1,^2) production is always associated with some amount 
of upstream momentum, either from initial state radiation or from decays of other, heavier 
particles. In this subsection we consider the effect of upstream momentum (Ut 7^ 0) and 
show that our previous conclusions still hold. 

It is well known that the kinematic endpoint M^°^ in general depends on the upstream 
momentum Ut, but this dependence is removed for a very special choice of the test masses, 
namely, when the test masses are equal to the true masses of the children particles, see 
eq. (4.7).^'^ This property has been previously demonstrated only for the case of identical 
parent particles and now we would like to test whether it also holds for the case of different 
parent particles in Fig. 6(a). 

As we already explained, the case of different parents in Fig. 6(a) can be equivalently 
treated as a case of identical parents, which decay asymmetrically as in Fig. 6(b). In turn, 
this process can be described in terms of the asymmetric Mt2D variable from eq. (4.1), 
whose endpoint M^'g(Mii,^ , Mi^/^) will become independent of Ut with the following choice 
of test masses 

Mvi>, = — ^ Mc, M^, = Mc. (5.18) 

These values are in principle measurable experimentally, by studying the Ut dependence 
of the M^2D endpoint as a function of the test masses, and finding the choice where this 
dependence is minimized [19]. 

Notice that if we try to use the symmetric version of the Mt2{Mc) variable, where the 
input masses are equal {Mci = = ^c), we will get a less stringent bound. Even if 
the test mass is taken to be the true one {Mq = Mq), we still find a chain of inequalities 

M^riMc) < M^iB Mc, Mc^ < Ma,. (5.19) 

Of course, when the test masses Mij,^ and Mi^,^ are chosen away from the special values 
in eq. (5.18), the endpoint M^2d ^ usual will be an increasing function of Ut- The same 
will also be true for its symmetric counterpart M^"^, whose Ut dependence is illustrated 



^"The dependence on Ut completely disappears for the case of the doubly projected Mt2± variable 
introduced in [40], which uses only the transverse momentum components orthogonal to Ut- 
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Figure 8. The same as Fig. 7, but for events with upstream momentum Ut- For the threshold 
hmit = + MA2) and for the infinite energy Hmit {V s — >■ 00), large Ut was put in by 
hand, while for the realistic LHC 8 TeV simulation, Ut was generated from initial state radiation 
in Pythia 6.4. 

in Fig. 8, which is the analogue of Fig. 7 in the presence of upstream momentum Ut- As 
before, the vertical blue and red dashed lines mark the locations of the expected endpoints 
in eqs. (5.16, 5.17) in the absence of Ut- The figure shows that, as expected, in the presence 
of upstream momentum, the endpoints are shifted higher, and in the extreme case of infinite 
^/s, the endpoint eventually reaches the naive expectation from eq. (5.1): 

lim M^2"'^(Mc = Mc) = Ma^- (5.20) 

Note that this limit is reached only in the unphysical case when both and Ut are sent to 
infinity. In the realistic simulation of squark-chargino production (the green dotted line), 
where Ut is generated from initial state radiation, the Mt2 endpoint is similar to the one 
observed in Fig. 7 and again does not saturate the bound eq. (5.20). 

5.3 Application to associated gluino-LSP production 

So far we have not at all discussed the event topologies in Figs.2(k) and 2(1), which can 
describe, e.g. associated gluino-neutralino production, where the gluino decay to the LSP 
gives 2 jets. Those topologies can be thought of as extreme examples of the "different 
parent" case just considered. Thus, M^"^' will again depend on the center of mass energy 
^/~^ of the AC2 system, which at hadron colliders will vary from one event to another. This 
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will lead to a smearing of the endpoint M-^"^. Thus we will consider the situation where 
the production energy \/5 of particles A and C2 is fixed, as at a future linear collider. 
The event topology of Fig. 2(k). The C/T-invariant point of M^'^^ will appear at 

(Mvp,, Mvt,) = {M^,, M^J , (5.21) 
M., = #ii-^fl-^V"l'^ (5.23) 



where 




rjb = cosh-i "^^^^ ^^^^ ^ ) . (5.25) 

The corresponding value of M^^^^^ will be 

M^'^^''\A'U„M^,) = ^. (5.26) 

The event topology of Fig. 2(1). With the effective topology technique, we find the Ut- 
invariant point at 

(^MvPi, Mvpj) = (M,i>, Mvt) , (5.27) 

with r]s still given by eq. (5.24). The endpoint is found at 

M^2D M^,) = MJp^'' (Mvp) = (5.29) 
6 The shapes of ^'1^2^ distributions 

Until now we have been focusing on the measurable kinematic endpoints of different vari- 
ables. At the same time, one could also attempt to study the shapes of the corresponding 
differential distributions. Unfortunately, (to the best of our knowledge) analytical formu- 
las for the shapes of the Mt2 and Mt2D distributions are absent. Their derivation would 
be rather complicated, because the shapes are affected by several factors: the production 
energy y/§, spin correlations, upstream momentum Ut, etc. In order to remove the Ut 
effect, Ref. [40] introduced a ID-projection of Mt2, called Mt2j^- Mt2^ is calculated the 
same way as Mt2, except that it uses the projections of the transverse momenta on the 
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dimension which is orthogonal to the Ut direction [49] . It turned out that this doubly 
projected variable is also independent of y/S and avoids large spin correlation effects [41]. 
In this section we provide analytical formulas for the shapes of the Mt2^ distributions for 
various event topologies from Fig. 2. We will also discuss the endpoint behavior of these 
distributions, extending the technique proposed in [26] to count the number of invisible 
particles Ninv- 

With two massless visible particles, Mt2^{M9) is related to Mj'2^(0) as follows 

The corresponding distributions are related as 

d / Mt2^ (m* ) 2 + M| \ d 



dMT2^(Afa.) I Mr2^(M„)2 / da; 



3.2) 



Thus we only need to describe the shape of x, since the shape of Mt2x(^^*) ^ill then be 
easily obtained from eq. (6.2). In the following equations, we used the same notation for ^ 
and rji as in eqs. (3.2, 3.4), respectively. We also introduce individual fii defined as 



Ma M? 
2 M2 



^^ = ^|1-t|^1> (6-3) 



so that fj, is the geometric mean of fii and fi2, as in eq. (3.2), 

M = VJmn- (6.4) 

In the following we list our results for dN/dx. 

The event topologies of Figs. 2(a), 2(b) and 2(e). The answer is very simple [40] 

dN /2/i\ , ^ 

— — oc X log — . (6.5) 

ax \ X J 

The event topologies of Figs. 2(c) and 2(f). We find 



AN dp 



(XX 

dx p 

2m 



41} (P) , (6.6) 



where 

4^! (P) = Vi-Q{p- e-^'O In {p e"' ) . (6.7) 

Here (x) is a unit step function and Jon {p) is a phase space weight from the cascade 
decay chain. After integrating out eq. (6.6), we get 



''2771 In f^^) if < X < 2ue-''i 



dA^ 

— — OC X • < 

dx 



if 2^e-''i < X < 2/i. 
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The event topologies of Figs. 2(d) and 2(g). We find 



dN 
dx 



OC X 



X q 

2p 



4^ il) 



\ 



{s-l)[s 



\Mc.-M^.\ 



(6.9) 
(6.10) 



The function J^g (q) is a three-body phase space weight. Integrating eq. (6.10) above, we 
get 



""^-"^^^-Ti^^yw- l^iog f M..v/^-M..7w^ 



Ma, \M^^^Ri-q^ + MA,^l-q\ 



log 



+ log V^,-g2 + ^l_g2 ^,(6.11) 



with 



When Xi is massless, 



Ma, = \Mc,-M^,\, Ri 



Ml - Ml 
Ml-Ml' 



(6.12) 



4^ i^) 



'Ml - Ml: ^ 



(l-g2)_ln 



Ml , 9s 9 



(6.13) 



and the corresponding distribution becomes 



diV 
dx 



OC X 



Ml-Mj.^ 



-l + -^-21n 



X 

2/x 



+Li2 



/2^ 

\Ma 



Li9 



2^2Ma 



(6.14) 



where Li2 (x) is Spence's function, defined as 

Li2 (x) = - 



(6.15) 



The event topology of Fig. 2(h). Without loss of generality, r/2 < rji and we have 



diV 
dx 



OC X • < 



In 

In 



if < X < xo. 



2/j 

xe 2 

m+V2 j 6r]iV2 
xe 2 



In 



2/j 



1 



r?2-31n(f +A 



In 



2/^ 



if Xo < X < xi, 



if Xi < X < X2, 



(6.16) 



In 



if X2 < X < X3, 
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Tchislepslep @LHC 8TeV, 80K events 



Endpoint: Topology (e) 




Mt2^(0) [GeV] 



Figure 9. The same as Fig. 5, but for the doubly projected Mt2 variable Mt2^- The solid lines 
show the corresponding theoretical predictions following from eqs. (6.5, 6.8, 6.16). 



where 



xo = 2/ie~('^^+''2\ xi = 2^ie-'^\ (6.17) 
X2 = 2iie-'^\ X3 = 2^i. (6.18) 



The event topology of Fig. 2(i). We find 



dN 

— — OC X • < 

dx 



if ? In (^) jS (<?) + /i- f (^i) j!S (g) , if < X < XI, 



(6.19) 



4?ln(^)j(J)(,), ifxi<X<X3. 
The event topology of Fig. 2(j). We leave it in integral form 

f „./:^,;<i.(„.S'(^). (0.0, 

For illustration, in Fig. 9 we show the corresponding Mt2± distributions for the Tchislepslep 
SMS model considered in Fig. 5. In addition to the histograms which were obtained from 
numerical simulations, we also show the corresponding theoretical predictions following 
from eqs. (6.5, 6.8, 6.16). We see that the analytical results match very well the simulated 
Mt2x differential distributions. 

As detailed in [26], it should, in principle, be possible to distinguish between topolo- 
gies with different numbers of invisible particles Ninv simply by fitting the endpoint fall-off 
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Table 1. Near-endpoint behavior of the doubly-projected stransverse mass Mt2± distribution, for 
all topologies in Fig. 2, as a function of e, the distance away from the endpoint. It is assumed 
that one invisible particle on each leg of the decay is massive (corresponding to the LSP), while 
any others are massless (corresponding to additional neutrinos emitted in the decay). The fall-off 
for alternative mass spectra can be trivially obtained by taking the relevant limit in the shape 
expressions given above. 

of kinematic distributions; this was shown to have a near- universal dependence on Ninv 
for judiciously chosen variables [26] . The near-endpoint behavior of the doubly-projected 
variable Mt2j_ for the various topologies in Fig. 2 is detailed in Table 1 as a function ofe, 
the distance away from the endpoint. Note firstly that the fall-off for the doubly-projected 
variable Mt2j^ is faster than that for the usual stransverse mass Mt2- This is due to an 
additional dimension of the full phase space being 'projected out', and making it more 
difficult to distinguish between different numbers of invisibles Ninv for large Ninv using the 
endpoint behavior alone. Secondly, although the near-endpoint behavior has a universal 
dependence on Ninv for massless invisible particles [26] , this is not true when any invisible 
particles obtain a non-negligible mass. In the latter case we see that the fall-off is always 
slower for true cascade decays (with on-shell resonances) , since the presence of an on-shell 
intermediate particle effectively reduces the dimension of the full phase space, hence fewer 
dimensions are projected out. In fact, for cascades, the near-endpoint behavior is entirely 
independent of the invisible particle masses. Finally, note that there are discrete ambigu- 
ities between the endpoint fall-off of cascades, and decays with off-shell intermediates, for 
different Ninv In these particular cases, then, one would need to fit using the full shape 
formula, given above. 

7 Conclusions and summary 

Our work in this paper removes some of the restrictions which so far have prevented the 
more widespread usage of the Mt2 variable. We demonstrated how the Mt2 variable (and 
its variants) can be usefully applied in more general situations, e.g.: 

• Decay chains with multiple invisible particles. In Sections 2 and 3 we considered cases 
where the new physics decay chain gives rise to several invisible particles. Previous 
M'T2 studies have typically assumed that there is only one invisible particle in the 
decay chain (the dark matter WIMP), which appears at the end of the decay chain. 
At the same time, there are many scenarios in which additional invisible particles 
can be present. The most popular example of this sort are chargino decays in super- 
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symmetry, which yield SM neutrinos in addition to the invisible LSP. We proposed 
two methods for dealing with the problem of additional invisibles: by introducing 
topology-dependent new variables (in Section 2) and by reinterpreting the measured 
kinematic endpoints of the conventional Mt2 distributions (in Section 3). 

• Events with different invisible child particles. In Section 4 we pointed out that the 
reinterpretation method carries over to the case where the invisible child particles 
at the end of the decay chains are different. The key idea is to use the asymmetric 
Mt2D variable introduced in [19]. 

• Events with different parent particles. Much of the previous literature on Mt2 dealt 
only with events in which the two parent particles initiating the decay chains are 
identical. As for events with different parent particles, it was thought that the M'T2 
endpoint in that case reveals the mass of the heavier parent. In Section 5 we showed 
that this conjecture is false, and we gave the correct interpretation of the Mt2 end- 
point for the case of different parent particles. 

Apart from a good theoretical understanding of the measured Mt2 endpoint in terms 
of the underlying mass spectrum in all those different situations, it is also important to 
have a good knowledge of the shapes of the respective differential Mt2 distributions. In 
Section 6 we considered the doubly projected variable Mt2j^ proposed in [40] and derived 
the corresponding shapes for a number of different cases shown in Fig. 2(a-j). Our formulas 
can be used for improving the precision of mass measurements based on Mt2^ kinematic 
endpoints [8]. Furthermore, by comparing the different shape predictions to the data, one 
could also, in principle, deduce the correct event topology and/or the number of invisible 
particles in the event, although in practice this may be unfeasible due to limited statistics 
near the endpoint. 
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